Stable fractional flux vortices and unconventional magnetic state in two-component 

superconductors. 
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In the framework of London theory we study the novel magnetic state in two-component super- 
conductors with finite density of fractional flux vortices stabilized near the surface. We show that 
the process of vortex entry into the two-component superconductor consists of several steps, while 
the external magnetic field increases from zero. At the first stage only vortices in one of the order 
parameter components penetrate and sit at the equilibrium position near the surface. When the 
magnetic field is increased further vortices in the second order parameter component eventually 
enter the superconductor. Such a complex partial vortex penetration leads to the modification of a 
Bean-Livingston barrier and a magnetization curve as compared to conventional single-component 
superconductors. We discuss the possibility of experimental identification of protonic superconduc- 
tivity in the projected superconducting state of liquid metallic hydrogen and hydrogen rich alloys 
with the help of partial vortex penetration effect. 

PACS numbers: 



I. INTRODUCTION 

Recently there has been a considerable interest 
to the physics of superconductors with multiple or- 
der parameter components, renewed by the discov- 
ery of a number of two-band superconducting mate- 
rials, such as MgB^, ferropnictidcs 2 -, heavy fcrmion 
compounds^, borocarbides^ and the advances in experi- 
mental techniques^ providing an intriguing possibility of 
observing new macroscopic quantum phenomena in liq- 
uid metallic hydrogen at ultra high pressures^—. Histori- 
cally the first realization of two component superconduct- 
ing state was considered by Moskalenko^ and by Suhl, 
Matthias, and Walker— in a metal with two overlapping 
energy bands on the Fermi surface. 

Multiple-component superconductors feature qualita- 
tively new effects, with respect to the conventional 
onesii. One of the striking differences is an unconven- 
tional structure of the mixed state in such superconduc- 
tors. For example, in Refill it was demonstrated that 
in two-band superconductor the interaction between vor- 
tices of the equal vorticity is not necessarily purely re- 
pulsive, but in some cases may be characterized by long- 
range attraction and short-range repulsion. This in par- 
ticular leads to the clusterization of vortices and forma- 
tion of the " semi-Meissner state", which was recently 
claimed to be experimentally observed in MgB^r^. 

Also a rich variety of peculiar vortex states in two- 
component superconductors have been investigated re- 
cently in the framework of the extended Ginzburg- 
Landau (GL) theory^. In particular, there was found 
a class of solutions of GL equations describing vortices 
bearing a fractional number of magnetic flux quanta. In 
general, such vortices can be characterized by different 
winding numbers of the order parameter in the two su- 
perconducting components, i.e. L\ ^ L2- Among such 
vortices one can distinguish a subclass of vortices which 
have L\ ^ and = 0, when the phase winding exists 



only in one of the condensates. In the present paper we 
will focus on the simplest case when L\ = ±1 and L-2 = 
(or equivalently L2 = ±1 and L\ =0), further using the 
term "fractional vortices" to identify such objects. Being 
positioned at the same point the two fractional vortices 
in different condensates characterized by the equal vor- 
ticity form a composite vortex bearing a single quantum 
of magnetic flux. 

The fractional vortices in two-component superconduc- 
tors are qualitatively different from thermodynamically 
stable exotic topological defects that occur in superfluid 
and superconducting systems with multiple order param- 
eter components, such us superfluid 3 He (see the review 
Ref. [H), ultracold atomsi& or vortices in d-wave HTSC 
cuprates with induced s-wave order parameter compo- 
nent inside vortex core^i. The fractional vortices consid- 
ered in the present work arc thermodynamically unstable 
in bulk two-component superconductor— since their en- 
ergy per unit length is logarithmically or linearly diver- 
gent with the sample size. Therefore it is impossible to 
create the fractional vortices in bulk superconductor by 
applying the external magnetic field. This is the reason 
why up to date fractional vortices in two component su- 
perconductors have not been observed. However it was 
proposed that the deconfinement of fractional vortices 
forming a composite vortex can occur due to the thermal 
fluctuation s 18 ' 19 or due to the thermal creation of frac- 
tional vortex-antivortex pairs^S. Also, with the help of 
GL calculations fractional vortices were demonstrated to 
exist in mesoscopic two-component superconductor s 21 ' 22 . 
In the present paper we propose that fractional vortices 
can be thermodynamically stable near the surface of su- 
perconductor and therefore there is a possibility to create 
a finite density of them by external magnetic field. The 
finite density of fractional vortices near the surface forms 
a novel magnetic state of superconductor which should 
be manifested by a modification of magnetization curve. 

Further we provide an analytical treatment of Lon- 
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don equations describing a behaviour of vortices near the 
surface of a two-component superconductor. We con- 
sider the two component superconductor as a mixture 
of two individually conserved superconducting conden- 
sates. This model is relevant for the systems where the 
interband coupling is forbidden by symmetry. Among 
such systems currently generating great interest is the 
liquid metallic hydrogen at ultra high pressures 5 -. In this 
case two superconducting components were predicted to 
originate from electronic and protonic Cooper pairing in 
metallic hydrogen and hydrogen rich alloys 6 - - — . Since the 
Cooper pairs of electrons cannot be converted into the 
Cooper pairs of protons the interband Josephson inter- 
action is strictly zero. An analogous model of two indi- 
vidually conserved condensates was considered recently 
in order to describe the exotic states of matter at neu- 
tron star inner cores with several charged barionic com- 
ponents, namely S _ hyperons and protons. The mixture 
of superfluid S - hyperons and protons is analogous to 
the two-gap superconductor with strictly zero interband 
Josephson couplin g 23 i 24 . 

On the other hand in two-band superconductors where 
the Cooper pairing of electrons takes place in differ- 
ent bands&^ the superconducting condensates in general 
cannot be considered as individually conserved. However 
the model with negligible interband Josephson coupling 
can also be applied to describe multiband superconduc- 
tors provided that all relevant physical scales are much 
shorter than the Josephson length. Although basically 
considering the system without the interband Josephson 
coupling we will discuss on a qualitative level the modi- 
fication of our results due to this effect. 

It is well known that the entry of vortices into a type-II 
superconductor is hindered by the so-called Bean- Liv- 
ingston surface barrie r 25 ' 26 . This barrier arises due to 
a competition of two forces acting on the vortex line: a 
force coming from the Meissner current driving the vor- 
tex into the superconductor and a force of the vortex 
mirror image attracting it towards the outside. As a re- 
sult, the penetration field H s of first vortex entry is typ- 
ically much larger than the lower critical field H c \ and 
becomes of the order of the thermodynamic critical field 
H<M. Recently it was demonstrated that a distribution 
of the Meissner current and the value of the penetration 
field H s are highly sensitive to the presence of Andrccv 
bound state at the surface of <i-wave superconductors 2 ^. 
It is the goal of the present paper to show that the two- 
component structure of the order parameter can also al- 
ter significantly the value of the penetration field H s and 
the process of the vortex entry in two-gap superconduc- 
tors. 

In two-component superconducting systems it is nat- 
ural to expect that the penetration of fractional vortices 
in the component with lower condensation energy should 
occur at magnetic field which is smaller than H c . How- 
ever, since the existence of fractional vortices is prohib- 
ited in the bulk, they should sit at a certain distance from 
the surface, corresponding to the energy minimum. The 



further increase of magnetic field will lead to the nuclc- 
ation of vortices in the second superconducting compo- 
nent. At a final stage of the process of vortex entry the 
fractional vortices of different types merge and proliferate 
into the bulk. 

This paper is organized as follows. In Sec. [H] we give 
an overview of the theoretical framework, namely the 
London theory of two-component superconductors. In 
Sec. IIIII we discuss the Gibbs energy of fractional vor- 
tices near the surface and address the questions of frac- 
tional vortex stability and the Bean-Livingston barrier 
for penetration of fractional vortices. Also we calculate 
the equilibrium distribution of fractional vortices near 
the surface. In Sec. II VI we discuss the modification of the 
magnetization curve due to the partial vortex penetra- 
tion and the possibility to implement the experimental 
identification of two order parameter components in the 
projected superconducting state of liquid metallic hydro- 
gen. We give our conclusions in Sec. [V] 



II. BASIC EQUATIONS: LONDON THEORY OF 
FRACTIONAL VORTICES IN TWO-GAP 
SUPERCONDUCTOR 

Let us consider a superconductor with two coexist- 
ing superconducting condensates. The external magnetic 
field is directed along the z axis Ho = HqZq. Consider- 
ing the London limit, i.e. assuming that the coherence 
length is vanishingly small compared to all other length 
scales we obtain the Gibbs energy per unit length along 
the z axis as follows^: 
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where the free energy is 
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Here Xa,b are the different length scales, proportional to 
the densities of two types of superconducting carriers, H 
is a total magnetic field, A is a vector potential, and 9 a,b 
are the phases of superconducting order parameters. As 
we will see further, the London limit in two-component 
model that we consider yields an exponential decaying 
magnetic field generated by fractional vortices. On the 
other hand it was recently demonstrated that in full two- 
component GL theory there is always a power-law tail 
of magnetic field around a fractional vortex^. However 
the power law behaviour of magnetic field starts at the 
distances greater than the London penetration depth A 
which is considered to be the largest length scale in our 
present paper. For example, as we will see below the 
intervortex distance in the proposed novel magnetic state 
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of two- component superconductor is much smaller than 
A. Within this parameter range the London model is a 
good approximation to study the magnetic properties of 
two-component superconductors. 

The London- Maxwell (LM) equation for the magnetic 
field obtained from the condition SF/SA = yields: 
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is the London penetra- 



</>o(A/Aa,s) 2 are partial mag- 



where A 2 = 
tion depth, and 4>a,b 
netic fluxes so that cf>A + 4>b = <fco- With the help of LM 
equation ([3]) the free energy (JT]) can be rewritten as a 
superposition of two parts 



F = F„ 



± nc 



(4) 



where the first part is the sum of total magnetic energy 
and the kinetic energy of superconducting current and 
thus depends only on the magnetic field 
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J [H 2 + A 2 (VxH) 2 ] d 2 r. 



(5) 



The second part F nc is the kinetic energy of relative mo- 
tion of two condensates, which can be associated with 
the neutral superfluid current, supported by the counter- 
directed motion of equally charged particles (electrons in 
two-gap superconductors) or codirected motion of oppo- 
sitely charged ones (electrons and protons in supercon- 
ducting liquid metallic hydrogen)ii. This term depends 
on the relative phase difference between two condensates 

frel = 9 A — Qb as follows 
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Such division of the total energy @ is convenient for the 
calculations since the parts F m and F nc can be found 
independently. 

For example the energy F m of a one fractional vor- 
tex in A or B condensate can be calculated exactly in a 
same scheme as for the case of conventional single gap 
superconductor—. The result is 



F m — Si + —H v (Ri), 

57T 



(7) 



where i = A, B and Ri is the fractional vortex coordi- 
nate. The first term in Eq.0 is the is the self energy of 
the vortex 

where & is the coherence length and the second term in 
Eq.© is the energy in the magnetic field H v = H v Zq 
generated by other vortices. 

At the infinite superconductng sample the part of the 
energy F nc given by ^ is divergent for a single frac- 
tional vortex due to the unscreened currents induced in 



the vortex free phased. Indeed, taking for example 9a = 
arctan(j//x) and 6b = we get that (V(^ re ;) 2 ~ 1/r 2 
which means that the expression ([6]) is logarithmically 
divergent with the size of the superconducting sample L 
so that F nc ~ InL. However further we will see that in 
some specific cases the energy of a single fractional vortex 
can be finite. In particular such situation is realized for 
fractional vortex placed near the boundary of supercon- 
ductor. In this case the large scale divergence in Eq.© is 
removed due to the cancellation of the unscreened super- 
conducting current in vortex free phase due to the image 
antivortex. 



III. RESULTS 

Gibbs energy of fractional vortices near the 
surface of superconductor. 



Image antivortices 



-XL 




A vortex B vortex 



-x 4 



FIG. 1: The sketch of the system under consideration. It 
consists of a two - component superconductor, occupying the 
half space x > bounded by the yz plane. The fractional 
vortices are situated at the distances xa and xb from the 
boundary. The boundary conditions are taken into account 
by placing image antivortices at the proper points behind the 
boundary plane. 

The sketch of the system considered is shown in Fig. 
[T] It consists of a two- component superconductor, oc- 
cupying the half space x > bounded by the yz plane. 
Our goal is to calculate the Gibbs energy ([T]) of the vor- 
tex configuration shown in Fig.l, i.e. of the two vortices 
in A and B phases located at the points Ra = (xa,0) 
and Rb = {xb, 0) near the surface. Here we can consider 
only vortices positioned in a line since any relative shift 
of A and B vortices along the boundary plane leads to 
the energy rising. 

All vortices are directed along the z axis, which is par- 
allel to the boundary plane. Since not Meissner currents 
nor vortices parallel to the surface do not create mag- 
netic fields outside the superconducotr— , the integration 
in Eq. ([T]) can be restricted to the superconducting region 
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only. Furthermore, in Eq.Q we can take into account 
only the magnetic field generated by vortices and neglect 
the field generated by the Meissner current. Then to cal- 
culate the energy we should find the magnetic field H(r) 
generated by vortices and the relative phase distribution 

^rei(r). 

The magnetic field is determined by the Maxwell- 
London equation having the following form 

A 2 VxVxH + H= (8) 
z [<f> A S{r - R A ) + <M(r - R B )] ■ 

Varying the expression (|6|) with respect to the phase dif- 
ference ip r ei we obtain the following 2D Poisson equation 
with the sources at the points of vortex positions 



A<p rel = 2ir [S(r - R A ) - S(r - R B )] 



(9) 



The Eqs. (|8l9|) should be supplemented by the bound- 
ary conditions at the surface of superconductor 
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\x=0 



0. 



(10) 



which means vanishing of the superconducting current 
through the boundary in both of the condensates. These 
boundary problems can be treated by using the method 
of images, i.e. by placing the image anti- vortices in A and 
B phases at the points Ra,b = (~xa,b,0) respectively 
[see Fig.©]. Then, the free energy of the two- vortex 
molecule near the flat surface of superconductor is just a 
one half of the free energy of vortex-antivortex molecule 
[Fig-©]. 

Let us evaluate the energy F m given by the expression 
([5]). For the stack of several vortices it is given by the 
sum of the individual vortex energies ([7]). The magnetic 
field H = H v zq generated by the vortex currents can be 
taken by the superposition of fields produced by vortices 
at the points Ra.b and antivortices at points Ra.b: 
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where Ko(x) is the zero - order Hankcl function, having 
the asymptotic K (x) w C — ln(x) at x <C 1 and Kq(x) 
v / n/2xe- x at x > 1 (RefH3). Thus for the energy F m 



we obtain: 
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where we introduce the coefficient a = cpA/tyB- 

To calculate the other part of free energy F nc we 
should substitute into Eq. © the relative phase distribu- 
tion in the form, corresponding to the case when A and 
B vortices are situated at the points Ra = (xa , 0) and 
Rb = (xbi 0) and image A and B antivortices are situ- 
ated at the points Ra = (— xa,0) and Rb = (— Xb,0) 
correspondingly: 
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|r-R A | 2 

The result is 
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Note that the expression (TTT|) for the energy F m as well 
as the expression (| 13[) for F nc are valid only when the sep- 
arations between vortices are larger than the coherence 
lengths: x A > £a, x b > and \x A ~~x B \ > max(^,^). 
However, the last restriction is removed when the total 
free energy 

(x A ,x B ) = F m (x A ,x B ) + F nc (x A ,x B ) 

is considered. Indeed, in this case the logarithmic singu- 
larities at xa = x B in Eqs. (TTT|) and (fT3")) cancel each 
other. Therefore the expression for the total free energy 
F vv {xa, x B ) given by the sum of partial energies (jlip and 
p3|) can be used for any values of the vortex coordinates 
x A > £a and x B > £ B . 

Finally let us consider the second term in Eq. (TTJ) which 
determines the interaction of vortices with the external 
magnetic field. This interaction energy can be written as 



W 



H $ 

Air ' 



(14) 



where $ = J H v d 2 r is the total magnetic flux gener- 
ated by the vortex currents. Due to the linearity of the 
London-Maxwell equation (JHJ and the boundary condi- 
tion (flQ|) the magnetic flux can be decomposed into the 
contributions from A and B vortices: 
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In accordance with the consideration of the Rcf . (1261) these 
partial magnetic fluxes depend on the positions of frac- 
tional vortices in the following way: 



the case of A phase vortex. Then the equation which 
determines x* A follows from the condition (f2"U)) and has 
the following form 



$ A = 



1 



cp B (l - e- XB ' x ) . 



(15) 



(16) 



Thus for the energy of vortex interaction with the exter- 
nal magnetic field we finally obtain: 

W(x A ,x B ) = 



Ho 

4-7T 



1 - er XA/x 



3 -xb/A 



(17) 



B. Stable fractional vortices 



Summing up the different parts of the free energy 
(|11I13P and the energy of interaction with the external 
magnetic field (fl7|) we obtain the Gibbs energy 



F G {x A ,x B ) = F vv (x Al x B ) + W(x A ,x B ). 



(18) 



The equilibrium vortex positions are determined by the 
cxtrcmum of the Gibbs energy, which is given by the 
condition 



OFg 
dxA 
dF G 
dx B 



(x A ,x* B ) — 
(x A , x* B ) = 



(19) 



if both x* A > £4 , x* B > £b or otherwise by the condition 



dF G 
dx A 
OFg 
dx B 



{x* A ,is) = Q 
{U,x* B ) = 



(20) 



when only one fractional vortex enters the superconduc- 
tor. 

As it is shown in left panel of Fig.@ for sufficiently 
large magnetic field there can be realized a situation when 
the fractional vortices are stable near the surface. In 
general it corresponds to the extremum condition (f2TJ)) 
which means that fractional vortices of different types can 
not coexist, i.e. either x* B — and x* A 7^ or x* A = and 
x* B 7^ 0. Another alternative is given by the condition 
(fT9")) which is satisfied only for a composite vortex with 
x* B = x* A 7^ 0, but it is impossible to realize a regime 
of stable molecule of fractional vortices near the surface 
with x* B x* A 7^ 0. 

Let us now discuss in more detail the dependence of the 
equilibrium position of a fractional vortex on the external 
magnetic field Hq = HqZq. To be definite we consider 



Ki(y) + — 



,y/2 
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where y = 2x* A /X and H£ = <j) A / \2tt\ 2 ). The solution 
of this equation is plotted in right panel Fig. J2} as a 
function of external magnetic field Hq for the different 
values of coefficient a = 1,2,3,4 (from top to bottom 
curve). Several qualitative features which characterize 
the behaviour of fractional vortices near the surface can 
be deduced from the Fig. ©b. Firstly, the distance x* A 
is scaled in the penetration depth A. This length scale 
is determined by the decay length of the interaction of 
vortex with the external field, i.e. the energy W(x A ). 
Secondly, the vortex coordinate x* A grows monotonically 
with the external field Hq. Such behaviour is explained 
by the fact that the force pushing the vortex towards the 
bulk of superconductor is proportional to Hq. Thus the 
larger is the pushing force the farther fractional vortex 
can penetrate into superconductor. Finally, the minimal 
magnetic field which provide the local Gibbs energy min- 
imum for x* A 7^ is of the order of — <p A / (2ttX 2 ). 

As we have shown above if the external magnetic 
field is larger than the threshold value of the order of 
H\ = 4>i/(2ir\ 2 ), where i = A,B there appears a min- 
imum of the fractional vortex energy which determines 
the equilibrium vortex position x* > 0. The Gibbs en- 
ergy minima occur both for the A and B vortices for an 
arbitrary ratio of the fractional fluxes <p A and 4> B as well 
as of the coherence lengths £ A and £g. For example in 
the left panel of Fig. @ the energy minima are shown to 
exist at x* A ^ 0, x* B = and x* B ^ 0, x* A = for a two 
gap superconductor with 4> A = 4> B and £4 = t; B . 

Thus we can conclude that fractional vortices can be 
stable near the surface of superconductor for arbitrary 
values of parameters <f> A , and <Pb,S,b characterizing the 
two superconducting condensates. Still the consideration 
above is not complete because it misses a very important 
question of how the fractional vortices can be created in 
the superconducting region. In the next section we will 
see that the scenario of vortex penetration is sensitive 
to the ratio of the coherence lengths £ A and t; B of two 
different condensates. 



C. Bean-Livingston barrier and the critical field of 
the first vortex entry 

To begin with the analysis of vortex penetration we 
note at first that the condition for a Bean-Livingston 
barrier suppression derived in Refl26l should be modified 
here to take into account the peculiarities of a two gap 
superconductor case. This condition takes the following 
form: 



(n- VF G ){x A =U,xb =£b) < 



(22) 
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different £4 7^ £b. In this case the vortex penetration 
occurs according to the foiiowing scenario. If the ex- 
ternal magnetic Hq field is increased from zero then at 
the threshold value Hq = H s i the Bean-Livingston bar- 
rier is suppressed and the fractional vortices with larger 
core size enter the superconductor. However, they can 
not proliferate into the bulk and sit at some equilibrium 
position near the surface of superconductor. The equi- 
librium position of one fractional vortex is determined 
by the Eq. (j2"Tj) and its dependence on the external field 
is shown in the right panel of FigJ^l When the mag- 
netic field Hq is increased further and reaches the second 
threshold field Hq = H S 2 determined by 



FIG. 2: Left panel: Gibbs energy of the molecule consisting 
of two fractional vortices near the surface of superconductor. 
The parameters are £a = £b and <f>A = 4>b- The magnetic 
field is taken Ho = 5H\. In this plot the local energy minima 
are clearly seen at is = 0, xa = ^a and xa = 0, xb = 
x* B . Right panel: The dependence of the fractional vortex 
position x* A on the magnetic field Ho for different values of 
4>a/4>b = 1,2,3,4 (from top to bottom curve). All lengths 
are normalized to the scale A. 



for some vector n = (n A ,n B ) satisfying the condition 
n A > 0, n B > 0. With the help of the expression for 
the Gibbs energy given by the Eqs.(HU [T3l H71 fT8]) it is 
possible to obtain an analytical expression for the critical 
value of magnetic field Hq = H s suppressing the surface 
barrier. 

It occurs that the two different regimes of vortex entry 
are determined by the ratio of coherence lengths of two 
superconducting components £4 and £b. 

(i) The case when £4 = £g = £ always corresponds 
to the simultaneous entrance of two fractional vortices. 
As one can see in Fig. ([3]) a the surface barriers for two 
vortices disappear at the same value of the critical field 
Hq = H s which is estimated as 



H s2 = max (H sAl H sB ) 



(25) 



H, 



4ttA£ 



(23) 



(ii) The case when £4 ^ £b is more interesting. In- 
deed in Fig.Qb it is clearly shown that there is a range 
of magnetic fields when for one of the superconducting 
condensates the Bean-Livingston barrier is already sup- 
pressed while for the other it still exists. The critical field 
of the first vortex entry Hq = H s i is given by 



H sl = mm(H sA ,H sB ) , 



(24) 



where H sA = (j>o/(AnX^ A ) and H sB = o /(47tA£b). 
Therefore the vortex with larger core size £a > £b is 
the first one to penetrate the superconductor. 

Thus we obtain that in some range of parameters it is 
possible that vortices penetrate by parts, i.e. only the 
fractional vortices in one of the condensates appear at 
first in the superconductor. Such situation is possible 
when the coherence lengths of the two condensates are 



the surface barrier is suppressed for the vortices in the 
second condensate. As one can see from Fig. ^ the 
fractional vortices in A and B condensates always merge 
to form a composite vortex and proliferate into the bulk 
superconductor since the total minimum of the Gibbs 
energy is always reached at x* A = x* B . Therefore we can 
conclude that if £4 7^ £b then at the range of magnetic 
fields H s i < Hq < H S 2 the fractional vortices should 
appear near the surface of a two-gap superconductor. 

Note that above we have analyzed only the single 
vortex problem. However, it is natural to expect that 
when the Bean-Livingston barrier for fractional vortices 
is suppressed they penetrate into the superconductor un- 
til some equilibrium vortex distribution is set up near 
the surface. Below we will find the equilibrium distribu- 
tion of fractional vortices at the range of magnetic fields 

H s i < Hq < H S 2- 



(b)4A=5 




FIG. 3: Gibbs energy of the molecule consisting of two frac- 
tional vortices near the surface of superconductor, (a) Equal 
coherence lengths of two condensates £a/£b = 1. In this 
case the Bean-Livingston barrier is suppressed simultaneously 
for fractional vortices of both types, (b) Different coherence 
lengths of two condensates £a/£,b =5. In this case the barrier 
is suppressed for the fractional vortex in A condensate while 
for the B condensate it still exists. 
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D. Unconventional magnetic state with finite 
density of fractional vortices near the surface. 

Now we proceed with considering the unconventional 
magnetic state of two-component superconductor which 
is realized when the finite density of stable fractional flux 
vortices appears near the surface as the magnetic field ex- 
ceeds the critical field of the first vortex entry H s i . To be 
definite we assume that £a > £,b and therefore the vortex 
in A phase is the first one to enter the superconductor 
since H sA < H sB . 

In general the task of finding the equilibrium config- 
uration of fractional vortices stabilized near the surface 
of superconductor seems to be a rather complicated one. 
One of the reasons is a many-body origin of this prob- 
lem, resulting in a huge number of dimensions of a con- 
figuration space. However we can use here a simplified 
approach which is based on the consideration of a distri- 
bution of an average density of fractional vortices. Such 
treatment is justified by the fact that in a magnetic field 
of the order H s a the intervortcx distance is of the order 
of a ~ y/X^A- As we have shown above the fractional 
vortices sit at the distance d ~ A from surface, which is 
much larger than the characteristic intervortex distance. 

So we assume that the A phase vortices are distributed 
in space with the density n(r), where r = (x,y). To find 
the form of the equilibrium vortex distribution n(r) we 
should consider a probe vortex positioned at some point 
r v = (x v , y v ) and place a requirement of its stationarity, 
i.e. that the total force acting on the probe vortex is to 
be zero 



d 

dx v 



F G {x v ) = 0. 



The expression for the Gibbs energy consists of three 
parts, namely 

F G (x v ) = F m {x v ) + F nc (x v ) + W(x v ), 

given by the Eqs. (JSJ |H [T4|) correspondingly. So to cal- 
culate the force we should determine the dependence of 
three energy parts on the position of the probe vortex. 
This leads to the following integral equation for the vor- 
tex density n(x) (the details of calculation are given in 
the Appendix |A"|) : 



(x X v ) ^_| z _ x<j | ^_\ X+Xv 
\x — x v \ 



n(x)dx+ (26) 



where a = ^> A /(16nX), (3 = §a<1>b /(87rA) and 7 



(3 / n(x)dx = 7e 



<PaF[o/(SttX). In Eq. (|2"6"|) the coordinates x and x v are 
normalized to the length scale A. 

The integral Eq. (|2f)|) can be solved analytically. Dif- 
ferentiating three times the both sides of Eq. (|2^|) over x v 
we obtain the differential equation for nix) 



dx 2 



13 + 2a 



(27) 



Therefore the solution n(x) has the form 

n(x) = Ae- x/L ", (28) 

where the scale L n and the amplitude A are determined 
by the Eqs. (|26l27l) as follows (here we restore the length 
normalization factor A): 



Ln = ^\ — 7T~ = AV 1 + 0a/4>b 



A = 1 (Ll 



A 2 ) _ H Q (Ll-X 2 ) 



2LI 



LI 



(29) 



(30) 



Thus we obtain that the concentration of fractional 
vortices n(x) given by Eqs. (|28l29l30[) decays exponen- 
tially with the distance from the surface of supercon- 
ductor. Note that in case of ordinary vortices (which is 
obtained by putting /3 = in the expressions above) we 
get that L n = 00, which means that the vortex density 
is constant^. The difference in the behaviour of ordi- 
nary and fractional vortex density n{x) is provided by 
an additional force which is determined by the gradient 
of energy F nc (x) ([6]). This force pushes the fractional 
vortices out of superconductor. Qualitatively this force 
is explained by the fact that the existence of fractional 
vortices is not allowed in the bulk superconductor due to 
the infinite energy of such objects. 

Let us find a contribution of fractional vortices to the 
average magnetization of a superconducting sample. To 
distinguish from the magnetization provided by Mcissncr 
currents we denote this contribution as Mf = ili/zo- 
Then we have: 



M 



1 



AttL 



n(x)$A(x)dx, 



(31) 



where n(x) is the vortex density, L is the sample size in 
x direction (see FigfTJ) and &a(x) is the total magnetic 
flux (fT5)) provided by one fractional vortex placed at the 
distance x from the boundary. Substituting expression 
(|28jl for the vortex density into Eq. t|3T|) we obtain that 



M f 



Ho L n — X 
4^r~ L ' 



(32) 



IV. DISCUSSION 



Now with the help of the results of previous section 
let us consider the properties of the magnetization curve 
of a two-gap superconductor. The qualitative features of 
the magnetization should be determined by the following 
facts. 

(i) In increasing magnetic field vortices penetrate into 
superconductor by parts. When the magnetic field 
reaches the critical value (f!R|) the fractional vortices with 
larger core size proliferate into superconducting region. 
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(ii) For the range of magnetic fields H s i < Hq < H s2 , 
only fractional vortices of one type exist near the surface 
of superconductor. Their distribution is determined by 
Eqs. (|28l29l30[) . When the magnetic field becomes larger 
Hq > H s2 the fractional vortices of the other type pen- 
etrate. Then fractional vortices of different types merge 
and composite vortices proliferate into the bulk supercon- 
ductor. The value of the field H s2 is certainly modified 
by the presence of finite density of fractional vortices in 
the first condensate. The modified value of H s2 can be 
found by calculating the forces forces acting of a single 
fractional vortex in the second condensate in a same way 
as on that in the first condensate (see AppendisTS} . Then 
we obtain th at in th is case the critical field is reduced by 
the factor y/4>B/<j>o as compared to Eq. ([25]) . where </>b 
is a magnetic flux carried by fractional vortices in the 
second condensate. 

Thus we can work out the following qualitative picture 
of the superconducting sample magnetization behaviour 
in increasing magnetic field. At low fields Hq < H s \ the 
magnetization is determined by the Meissner current. If 
all the stray fields are neglected it is given by 



Ah 



Ho 

4-7T ' 



(33) 



For magnetic fields from the range H sl < Hq < H s2 
the magnetization is changed due to the fractional vor- 
tices 



M 2 = M f (H ) - 



(34) 



where Mf(H ) is given by Eq.(|3"2]). Finally at H > H s2 
the magnetization changes abruptly and becomes 



M 3 = 



B(Hq) - Hq 

4tt 



(35) 



where the net magnetic induction B(Hq) is determined 
by the configuration of the vortex lattice in bulk super- 
conductor (see for example Refs l26lf29| ). 

These qualitative features of the magnetization curve 
are summarized in Fig.(j4]), where the magnetic field de- 
pendence of the average magnetization of superconduct- 
ing sample M(Hq) is shown for the case of single gap 
(left panel) and two gap (right panel) superconductors. 
For a single-gap superconductor we have a conventional 
picture of the Meissner state overheating due to the 
Bean-Livingston barrier. In the left panel of Fig. ((4]) we 
show the three characteristic regions corresponding to 
this case. At the region (I) when the magnetic field is 
smaller than the first critical one H c \ the Meissner state 
is realized. At the region (II) where H c \ < Hq < H s , 
where H s = 0o/(47rA£) the Meissner state is overheated 
due to the Bean-Livingston barrier. Correspondingly, the 
total magnetization still grows at this region (solid line) 
instead of following the bulk sample curve (dash line). 
Finally at Hq = H s the surface barrier is suppressed and 
the magnetization jumps to the value determined by the 
configuration of vortices in the bulk of superconductor. 



In contrast to the conventional scheme a two gap super- 
conductor features two jumps of the magnetization curve 
[right panel of Fig.(j4|)]. Correspondingly there are four 
regions which are characterized by qualitatively different 
behaviour of the M{Hq) dependence. The regions (I) 
where Hq < H c \ and (II) where H c \ < Hq < H s i arc the 
same as for conventional superconductor. A nontrivial 
behaviour starts at the field Hq = H s \ of the first vortex 
entry. At this threshold field value there is a jump of the 
magnetization curve, which is determined by a set up of 
a finite concentration of fractional vortices near the sur- 
face of superconductor. At the region (III) for the range 
of magnetic fields H s i < H < H s2 the magnetization is 
determined by the Eqs. (|32I34I) and grows linearly with 
Hq. The second jump of the M(Hq) curve occurs at the 
field Hq = H s2 when the composite vortices start to pro- 
liferate into the bulk of superconductor. Finally at the 
region (IV) when Hq > H s2 the value of magnetization 
is determined by the configuration of composite vortices 
in the bulk of superconductor. 

The discussed above possibility of separate vortex pen- 
etration in two coexisting superconducting condensates 
should be especially interesting in connection with the 
investigation of projected superconducting state of liquid 
metallic hydrogen^, where the superconducting state is 
formed by electronic and protonic Cooper pairs. The ob- 
servation of the two-component superconducting state in 
this case cannot be implemented by conventional tech- 
niques and requires special experimental approaches^. 
In particular, the most challenging problem of the pro- 
tonic superconductivity detection cannot be treated by 
the standard measurement of the Meissner effect since the 
critical temperature for electrons is estimated to be much 
larger than that of protons. Therefore, in the Meissner 
state the contribution to the total magnetic moment of 
protonic supercurrent will be always masked by that of 
the electronic component. On the other hand, as we have 
shown above the relaxation of the overheated Meissner 
state by vortex penetration should feature an additional 
jump of magnetization due to the coexistence of protonic 
and electronic superconducting components. According 
to our results the critical magnetic field of the first vortex 
entry is determined by the condensate with the largest 
coherence length, which is the one with the smallest crit- 
ical temperature. In liquid metallic hydrogen such com- 
ponent is always a protonic one, therefore the first jump 
on the magnetization curve should be determined by the 
vortices in protonic supcrfluid provided it is a type-II one. 
According to the recent estimations the low temperature 
limit of magnetic flux carried by vortices in protonic com- 
ponent is of the order (f> p (T = 0) ~ 10~ 3 </>o being by 
several orders of magnitude larger than the resolution 
threshold in modern experiments^ which allow to detect 
the magnetic flux less then lO _5 0o- Therefore even for 
the temperatures close to the critical one of protinic com- 
ponent it is possible to detect the protonic vortices. It is 
interesting to note that on approaching the critical tem- 
perature of protonic component from below the penetra- 
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tion field for the fractional vortex in protonic condensate 
goes to zero. Therefore the largest interval of magnetic 
fields where only the fractional vortices should exist near 
the surface is in the vicinity of protonic critical temper- 
ature. 

Finally we note that besides an intriguing possibility 
to explore the nature of superconducting state in liquid 
metallic hydrogen our results are applicable to the case of 
conventional two-gap superconductor. Generally speak- 
ing in this case one should take into account the effect 
of interband Josephson coupling that we have neglected. 
This assumption is justified in case when the Josephson 
length of the intercomponent phase difference relaxation 
is much larger than the other relevant length scales, that 
are the coherence lengths in both condensates and the 
London penetration length. However even if this con- 
dition is violated, the critical magnetic field of the first 
vortex entry H s \ remains the same, i.e. it is not af- 
fected by the presence of the Josephson coupling since 
the energy of Josephson string connecting the fractional 
vortex with the surface grows quadratically at small dis- 
tances (see e.g. Refl3ll) and does not alter the condition 
of vortex penetration (|22[) . The further penetration of 
vortex in the second component is certainly affected by 
the presence of the first fractional vortex emanating the 
Josephson string ending at the surface. Presumably, the 
result should be the reduction of the field H S 2 of the sec- 
ond vortex entry, however the strict quantitative inves- 
tigation of this effect is beyond the scope of the present 
paper. Thus we can conclude that our main result, i.e. 
the separate penetration of fractional vortices and the 
presence of two jumps on the magnetization curve re- 
mains qualitatively relevant even in case of rather strong 
intercomponent Josephson interaction. 
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FIG. 4: Qualitative behaviour of magnetization curve of a su- 
perconductor modified due to the surface barrier. Left panel: 
Single gap superconductor. Right panel: Two-gap supercon- 
ductor. 



and therefore can not exists. In our situation the sta- 
bility of fractional vortices near the surface is provided 
by the cancellation of the unscreened superconducting 
current due to the image antivortices. 

Also we have discussed the influence of fractional vor- 
tices on a Bean-Livingston barrier and the magnetization 
curve of a two-component superconductor. In particu- 
lar, we have found that if the coherence lengths of two 
condensates are different then vortices penetrate the su- 
perconductor by parts. The fractional vortices in the 
condensate with larger coherence length, i.e. having the 
larger core size are the first to enter the superconductor. 
By increasing the magnetic field further the fractional 
vortices in the second condensate are pushed to pene- 
trate. Then the fractional vortices of two different types 
merge to form composite vortices which proliferate into 
the bulk superconductor. 

We have shown that the magnetization curve of a su- 
perconduting sample should feature the two jumps as- 
sociated with the penetration of two types of fractional 
vortices. The observation of such a peculiar magnetiza- 
tion behaviour could be considered as an experimental 
identification of fractional vortices in systems with sev- 
eral superconducting components, such as multiband su- 
perconductors and the projected superconducting state 
of liquid metallic hydrogen. 
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Appendix A: Derivation of the integral equation (|2(i[) 



V. SUMMARY 

To sum up we have investigated the unconventional 
magnetic state of two-component superconductors. This 
peculiar state is realized when the finite density of stable 
fractional vortices appears near the surface of supercon- 
ductor under the action of external magnetic field. This 
result contrasts to the case of bulk two-component super- 
conductor where fractional vortices have infinite energy 



(i) To find the energy F m (x v ) let us consider at first 
the interaction of the probe vortex with vortex positioned 
at the coordinate r = (x, y) and the antivortex at the 
coordinate r = (—x, y). Applying the Eq.© we get into 
account only the energy of interaction between vortices 
and obtain 

F m (r v ,r) = (Al) 
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Now to find the total energy F m (x v ) we integrate over 
the positions of all vortices to obtain: 



The integration over the coordinate y in Eq. (|A5|) can 
be performed using the relation 



-^m {Xy ) 



dy / dx n(x)F m (r v ,r). 



(A2) 



To perform the integration in Eq. (|A2j) we use the follow- 
ing relation 



Kq (yjx 2 + y 2 ^j dy = ire |a:| 



and obtain 

{Xv ) 



8A?r J 



exp 
exp 



\X — Xy\ 

A 

\x + x v \ 
A 



(A3) 



n(x)dx. 



(ii) According to the Eq. © energy F nc (x v ) is deter- 
mined by the following expression 



Fnc {.Xy ) 



VA<PB 



tt (4ttA) 2 



(A4) 



where 



= arctan 



V - Vv 



is the phase distribution created by the probe vortex and 

y-Vi 



E 



to, arctan 



X Xi 



is the phase created by all other vortices (to^ = 1) and 
antivortices (to.; = — 1). To evaluate the expression (|A4[) 
let us consider again the interaction of the probe vortex 
with vortex positioned at the coordinate r = (x, y) and 
the antivortex at the coordinate r = (—x,y). Then we 
have that 



F nc (r V! r) = 2^1n / 



(4ttA) 2 



r - r v 



This expression is divergent at the point r = r v , which 
should be cut off at the vortex core size However, 
this singularity is integrablc and we don't introduce this 
cutoff here. Then we should sum the contributions from 
all vortices according to the expression 



dy I dx n(x)F nc (r v , r). (A5) 



||r + r v |\ _ (2nx v for x > x v 
|r — r„| / 1 2ttx for x < x„ 



Then we obtain: 

F nc (x v ) 

4>a4>b 

4ttA 2 



(A6) 



xn(x)dx + / x v n{x)dx 



(iii) Finally the energy W(x v ) determining the inter- 
action of the probe vortex with external magnetic field 
can be find from Eqs. (|14ll5j) : 



W(x v ) = -^H (: 



(A7) 



Summing all the contributions to the Gibbs energy 
()A3IA6|A7j) we obtain the total force acting on a probe 
vortex: 



d 

f = g^r^G = (/to + fnc + ./ezt) X 0, 



where 



fm. {Xv ) 

+ exp 



\x + x v \ 
A 



(x - x v ) 
\x - x v \ 

n(x)dx, 



exp 



\x + x v \ 
A 



(A8) 



and 



fnc = / n(x)dx 



f _ <PA J] -X V /X 



(A9) 



(A10) 



The condition of the probe vortex stationarity f = 
yields the integral Eq. ((26|) . where the coordinates x and 
x v are normalized to the length scale A. 
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